Sensitivity analysis and constraint of the parameter space of Earth System configuration of JULES

1 Introduction

Update of explore-JULES-ES-1p0.Rmd that uses the new packages and functions

We build and test gaussian process emulators, perform a sensitivity analyis, and constrain the input space of Jules.

Andy thinks there might be mileage in analysing the atmospheric growth, which is here. /home/h01/hadaw/Research/ES_PPE/Oct20

At the moment, this vignette is hampered by the fact that emulators are failing on a few of the outputs which represent change over the historical period. The emulator is fine for predicting absolute values in the modern period.

Andy would like to see timeseries of: cVeg, cSoil and nbp, npp in GtC and GtC/yr.

1.1 Preliminaries

Load libraries, functions and data.

# Load helper functions

knitr::opts_chunk$set(fig.path = "figs/", echo = FALSE, message = FALSE, warnings = FALSE)
# load some helper functions
source('~/brazilCSSP/code/brazil_cssp/per_pft.R') # eventually, move the relevant functions
source('explore-JULES-ES-1p0_PPE_functions.R')

1.2 Carbon budget data

Section 2.5 in Friedlingstein et al. describes how the land carbon sink is estimated.

1.2.1 Jules output data location

1.3 Extract a time series of an annual, globally aggregated variable

1.4 Plot timeseries of aggregated global timeseries

1.5 Can we understand change and absolute?

1.6 Function to extract the “modern value” direct from the file (last 20 years of the timeseries)

1.7 Loop to extract the “modern value” of a number of model outputs

1.8 Calculate an ensemble of anomalies for all variables

For each ensemble member and each variable, calculate the change from the 20 years at the start of the run, to the twenty years at the end of the run.

1.9 Clean data sets to “level 0”

Initial clean of data set, removing variables that don’t change, and removing NAs (models that didn’t run).

[1] TRUE

1.10 Where did the model fail to run?

(“probability of run failure” Could be an interesting project - logit transformation for probability of run failure? Some other ML like SVM?

There are clear run failure thresholds in the parameters rootd_ft_io and lai_max, and quite strong visual indications that a_wl_io and bio_hum_cn matter.

NEXT, histograms of failure with the LHS value

Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory

Error in nc_open(paste0(fn)) : 
  Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc

Error in R_nc4_inq_varndims: NetCDF: Not a valid ID

Error in ncvar_ndims(ncid, varid) : error returned from C call

Error in R_nc4_close: NetCDF: Not a valid ID

1.11 Testing Gaussian Process emulators

1.11.1 NPP as a test

First, use mean NPP as an example. How does NPP respond to each parameter? NAs are removed, but zero values are still included.

1.12 Relationship between modern NPP and change since 1850.

Some outputs (e.g fLuc, fHarvest) have an almost perfect 1:1 relationship between modern value and change, some (nbp, npp, treeFrac) quite or moderately strong, and some (csoil, cveg) very weak or non-existant.

1.13 A clear threshold in the F0 parameter for NPP

It appears that this ensemble is less “clear cut” in having an output that clearly distinguishes between “failed” (or close to it), and “not failed”.

Having said that, having an F0 over a threshold seems to kill the carbon cycle, as before. Here, we’ve set a threshold of 0.9 (on the normalised scale) for F0, and we remove members of the ensemble with a larger F0 than that when we build emulators.

1.14 Level 1 constraint - remove ensemble members with F0 above threshold

The level 1 constraint removes any input with F0 greater than 0.9 (normalised), which removes many of the zero-carbon-cycle members up front. There are 424 ensmble members remaining.

This does constraint sequentially, which may not be a good idea.

[1] 2.410932e-09

1.14.1 Comparison of level 0 and level 1 emulators for NPP (modern value)

# Plot the regular km emulator. Doesn’t look great.

1.14.2 Leave-one-out summaries of NPP emulators

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2 
  - best initial criterion value(s) :  -1842.508 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1842.5  |proj g|=       1.8652
At iterate     1  f =       1841.3  |proj g|=        1.4365
At iterate     2  f =       1839.3  |proj g|=        1.8697
At iterate     3  f =         1839  |proj g|=        1.8947
At iterate     4  f =         1839  |proj g|=        1.2143
At iterate     5  f =       1838.6  |proj g|=        1.1891
At iterate     6  f =       1838.3  |proj g|=        1.9015
At iterate     7  f =       1837.6  |proj g|=         1.871
At iterate     8  f =         1837  |proj g|=        1.8711
At iterate     9  f =       1834.9  |proj g|=        1.7796
At iterate    10  f =         1833  |proj g|=        1.8503
At iterate    11  f =       1831.1  |proj g|=        1.8644
At iterate    12  f =       1830.9  |proj g|=        1.8531
At iterate    13  f =       1829.9  |proj g|=        1.8044
At iterate    14  f =       1828.8  |proj g|=        1.8322
At iterate    15  f =       1827.7  |proj g|=        1.8111
At iterate    16  f =       1825.8  |proj g|=        1.8566
At iterate    17  f =       1822.1  |proj g|=        1.8783
At iterate    18  f =       1817.9  |proj g|=        1.8192
At iterate    19  f =         1810  |proj g|=        1.8972
At iterate    20  f =       1807.7  |proj g|=        1.3523
At iterate    21  f =       1806.7  |proj g|=        1.3328
At iterate    22  f =       1804.9  |proj g|=        1.8829
At iterate    23  f =       1803.2  |proj g|=        1.6989
At iterate    24  f =       1802.4  |proj g|=        1.7693
At iterate    25  f =       1801.2  |proj g|=        1.8451
At iterate    26  f =       1799.8  |proj g|=        1.8602
At iterate    27  f =       1799.5  |proj g|=        1.6174
At iterate    28  f =       1799.2  |proj g|=        1.5947
At iterate    29  f =       1798.8  |proj g|=        1.5503
At iterate    30  f =       1796.6  |proj g|=        1.8546
At iterate    31  f =       1795.7  |proj g|=         1.586
At iterate    32  f =       1795.4  |proj g|=         1.583
At iterate    33  f =       1795.3  |proj g|=        1.6047
At iterate    34  f =       1793.7  |proj g|=        1.5063
At iterate    35  f =       1793.4  |proj g|=        1.4599
At iterate    36  f =       1793.2  |proj g|=        1.8832
At iterate    37  f =         1793  |proj g|=          1.87
At iterate    38  f =       1792.9  |proj g|=        1.4333
At iterate    39  f =       1792.4  |proj g|=        1.6437
At iterate    40  f =       1792.1  |proj g|=        1.3598
At iterate    41  f =       1791.8  |proj g|=        1.8392
At iterate    42  f =       1791.2  |proj g|=        1.8717
At iterate    43  f =       1790.7  |proj g|=        1.8739
At iterate    44  f =       1790.5  |proj g|=        1.0022
At iterate    45  f =       1790.4  |proj g|=        1.1185
At iterate    46  f =       1790.3  |proj g|=        1.8759
At iterate    47  f =       1790.3  |proj g|=         1.875
At iterate    48  f =       1790.2  |proj g|=        1.8727
At iterate    49  f =       1790.1  |proj g|=         1.802
At iterate    50  f =       1790.1  |proj g|=       0.81547
At iterate    51  f =       1790.1  |proj g|=        0.8373
At iterate    52  f =         1790  |proj g|=        1.8728
At iterate    53  f =       1789.8  |proj g|=       0.93276
At iterate    54  f =       1789.8  |proj g|=        1.8773
At iterate    55  f =       1789.7  |proj g|=        1.8754
At iterate    56  f =       1789.7  |proj g|=        1.5203
At iterate    57  f =       1789.6  |proj g|=        1.2658
At iterate    58  f =       1789.6  |proj g|=       0.92604
At iterate    59  f =       1789.2  |proj g|=         1.876
At iterate    60  f =       1789.2  |proj g|=         1.873
At iterate    61  f =       1789.1  |proj g|=       0.80926
At iterate    62  f =       1789.1  |proj g|=       0.63617
At iterate    63  f =       1789.1  |proj g|=        1.7888
At iterate    64  f =         1789  |proj g|=         1.873
At iterate    65  f =         1789  |proj g|=        1.4746
At iterate    66  f =         1789  |proj g|=       0.71729
At iterate    67  f =         1789  |proj g|=       0.71538
At iterate    68  f =         1789  |proj g|=       0.70747
At iterate    69  f =         1789  |proj g|=       0.63382
At iterate    70  f =       1788.9  |proj g|=       0.74475
At iterate    71  f =       1788.9  |proj g|=        1.7914
At iterate    72  f =       1788.9  |proj g|=       0.95343
At iterate    73  f =       1788.8  |proj g|=       0.85278
At iterate    74  f =       1788.8  |proj g|=       0.74093
At iterate    75  f =       1788.8  |proj g|=       0.90557
At iterate    76  f =       1788.8  |proj g|=       0.29949
At iterate    77  f =       1788.8  |proj g|=       0.28806
At iterate    78  f =       1788.8  |proj g|=       0.27628
At iterate    79  f =       1788.8  |proj g|=       0.55663
At iterate    80  f =       1788.8  |proj g|=        1.7132
At iterate    81  f =       1788.8  |proj g|=        1.8713
At iterate    82  f =       1788.8  |proj g|=        1.8716
At iterate    83  f =       1788.8  |proj g|=        1.8729
At iterate    84  f =       1788.8  |proj g|=       0.83666
At iterate    85  f =       1788.8  |proj g|=       0.30137
At iterate    86  f =       1788.7  |proj g|=       0.31882
At iterate    87  f =       1788.7  |proj g|=        0.2255
At iterate    88  f =       1788.7  |proj g|=       0.42148
At iterate    89  f =       1788.7  |proj g|=       0.59752
At iterate    90  f =       1788.7  |proj g|=       0.51451
At iterate    91  f =       1788.7  |proj g|=        1.4762
At iterate    92  f =       1788.7  |proj g|=       0.69297
At iterate    93  f =       1788.7  |proj g|=       0.53596
At iterate    94  f =       1788.7  |proj g|=        0.5359
At iterate    95  f =       1788.7  |proj g|=       0.55123
At iterate    96  f =       1788.7  |proj g|=       0.37871
At iterate    97  f =       1788.7  |proj g|=       0.31548
At iterate    98  f =       1788.7  |proj g|=       0.23863
At iterate    99  f =       1788.7  |proj g|=       0.13991
At iterate   100  f =       1788.7  |proj g|=       0.13728
At iterate   101  f =       1788.7  |proj g|=       0.14306
final  value 1788.683595 
stopped after 101 iterations

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2 
  - best initial criterion value(s) :  -6192.077 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       6192.1  |proj g|=       1.7779
At iterate     1  f =       6097.2  |proj g|=        1.9547
At iterate     2  f =       6092.5  |proj g|=        1.9408
At iterate     3  f =       6091.7  |proj g|=          1.94
At iterate     4  f =       6089.5  |proj g|=        1.9104
At iterate     5  f =       6081.3  |proj g|=        1.9157
At iterate     6  f =       6077.6  |proj g|=        1.9116
At iterate     7  f =       6074.3  |proj g|=         1.907
At iterate     8  f =       6073.6  |proj g|=        1.9049
At iterate     9  f =       6071.5  |proj g|=        1.8255
At iterate    10  f =         6066  |proj g|=        1.7837
At iterate    11  f =       6065.2  |proj g|=        1.7673
At iterate    12  f =       6063.4  |proj g|=         1.752
At iterate    13  f =       6062.6  |proj g|=        1.9208
At iterate    14  f =       6062.1  |proj g|=        1.7448
At iterate    15  f =         6062  |proj g|=        1.7437
At iterate    16  f =       6061.8  |proj g|=        1.7402
At iterate    17  f =         6061  |proj g|=        1.7269
At iterate    18  f =       6060.1  |proj g|=        1.7134
ys=-1.254e-02  -gs= 8.259e-01, BFGS update SKIPPED
At iterate    19  f =       6059.6  |proj g|=        1.7043
ys=-2.746e-03  -gs= 5.566e-01, BFGS update SKIPPED
At iterate    20  f =       6059.3  |proj g|=        1.7005
At iterate    21  f =       6056.3  |proj g|=        1.9169
At iterate    22  f =       6055.7  |proj g|=        1.9087
At iterate    23  f =       6055.4  |proj g|=        1.8996
At iterate    24  f =       6055.1  |proj g|=        1.6244
At iterate    25  f =       6054.9  |proj g|=        1.8958
At iterate    26  f =       6054.6  |proj g|=        1.6078
At iterate    27  f =       6054.3  |proj g|=        1.5928
At iterate    28  f =       6053.9  |proj g|=        1.8949
At iterate    29  f =       6053.8  |proj g|=        1.8942
At iterate    30  f =       6053.7  |proj g|=        1.5209
At iterate    31  f =       6053.7  |proj g|=         1.516
At iterate    32  f =       6053.3  |proj g|=        1.8931
At iterate    33  f =       6052.9  |proj g|=        1.3806
At iterate    34  f =       6052.8  |proj g|=        1.3433
At iterate    35  f =         6052  |proj g|=        1.8927
At iterate    36  f =       6051.9  |proj g|=        1.3527
At iterate    37  f =       6051.9  |proj g|=        1.0531
At iterate    38  f =       6051.8  |proj g|=       0.96661
At iterate    39  f =       6051.8  |proj g|=       0.97303
At iterate    40  f =       6051.8  |proj g|=       0.62567
At iterate    41  f =       6051.7  |proj g|=       0.68667
At iterate    42  f =       6051.7  |proj g|=        1.0212
At iterate    43  f =       6051.7  |proj g|=       0.36663
At iterate    44  f =       6051.7  |proj g|=       0.20017
At iterate    45  f =       6051.7  |proj g|=       0.20136
At iterate    46  f =       6051.7  |proj g|=        0.1629
At iterate    47  f =       6051.7  |proj g|=       0.21978
At iterate    48  f =       6051.7  |proj g|=       0.23363
At iterate    49  f =       6051.7  |proj g|=       0.17448
At iterate    50  f =       6051.7  |proj g|=      0.091728
At iterate    51  f =       6051.7  |proj g|=      0.097416
At iterate    52  f =       6051.7  |proj g|=       0.12446
At iterate    53  f =       6051.7  |proj g|=       0.07705
At iterate    54  f =       6051.7  |proj g|=       0.11677
At iterate    55  f =       6051.7  |proj g|=       0.20318
At iterate    56  f =       6051.7  |proj g|=      0.094594
At iterate    57  f =       6051.7  |proj g|=       0.11668
At iterate    58  f =       6051.7  |proj g|=       0.03121
At iterate    59  f =       6051.7  |proj g|=      0.055344
At iterate    60  f =       6051.7  |proj g|=       0.11291
At iterate    61  f =       6051.7  |proj g|=      0.066099
At iterate    62  f =       6051.7  |proj g|=      0.078269

iterations 62
function evaluations 76
segments explored during Cauchy searches 76
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0782694
final function value 6051.72

F = 6051.72
final  value 6051.720344 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2 
  - best initial criterion value(s) :  -1844.124 -1844.169 -1844.496 -1844.519 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1844.1  |proj g|=       1.9783
At iterate     1  f =       1827.3  |proj g|=        1.5081
At iterate     2  f =       1810.4  |proj g|=        1.9257
At iterate     3  f =       1810.1  |proj g|=        1.2762
At iterate     4  f =       1808.6  |proj g|=        1.9155
At iterate     5  f =       1803.8  |proj g|=        1.9384
At iterate     6  f =       1802.5  |proj g|=        1.9352
At iterate     7  f =       1801.9  |proj g|=        1.9195
At iterate     8  f =       1801.3  |proj g|=        1.5673
At iterate     9  f =         1801  |proj g|=        1.5746
At iterate    10  f =       1800.3  |proj g|=        1.5734
At iterate    11  f =         1800  |proj g|=        1.5698
At iterate    12  f =       1799.2  |proj g|=        1.9053
At iterate    13  f =         1799  |proj g|=        1.9022
At iterate    14  f =       1798.9  |proj g|=        1.2078
At iterate    15  f =       1798.9  |proj g|=        1.2015
At iterate    16  f =       1798.8  |proj g|=        1.4287
At iterate    17  f =       1798.3  |proj g|=        1.9001
At iterate    18  f =       1797.6  |proj g|=        1.8921
At iterate    19  f =       1796.9  |proj g|=        1.8927
At iterate    20  f =       1796.5  |proj g|=        0.9502
At iterate    21  f =       1795.9  |proj g|=        1.2526
At iterate    22  f =       1795.6  |proj g|=        1.8954
At iterate    23  f =       1795.4  |proj g|=        1.5256
At iterate    24  f =       1794.9  |proj g|=        1.8901
At iterate    25  f =       1794.7  |proj g|=        1.5025
At iterate    26  f =       1793.9  |proj g|=        1.3057
At iterate    27  f =       1793.7  |proj g|=        1.3204
At iterate    28  f =       1793.5  |proj g|=        1.4504
At iterate    29  f =       1793.4  |proj g|=        1.2978
At iterate    30  f =       1793.3  |proj g|=        1.2918
At iterate    31  f =       1793.1  |proj g|=        1.8873
At iterate    32  f =       1792.8  |proj g|=        1.8811
At iterate    33  f =       1792.6  |proj g|=        1.8802
At iterate    34  f =       1792.2  |proj g|=        1.7755
At iterate    35  f =       1791.9  |proj g|=        1.7764
At iterate    36  f =       1790.9  |proj g|=        1.7796
At iterate    37  f =       1790.8  |proj g|=        1.8808
At iterate    38  f =       1790.7  |proj g|=        1.2525
At iterate    39  f =       1790.7  |proj g|=       0.78214
At iterate    40  f =       1790.4  |proj g|=       0.82047
At iterate    41  f =       1790.3  |proj g|=        1.0794
At iterate    42  f =       1790.2  |proj g|=        1.3994
At iterate    43  f =       1790.1  |proj g|=        1.8787
At iterate    44  f =       1789.9  |proj g|=        1.2268
At iterate    45  f =       1789.9  |proj g|=       0.90265
At iterate    46  f =       1789.8  |proj g|=        1.8777
At iterate    47  f =       1789.8  |proj g|=        1.8783
At iterate    48  f =       1789.7  |proj g|=        1.5976
At iterate    49  f =       1789.7  |proj g|=       0.96825
At iterate    50  f =       1789.7  |proj g|=        1.7855
At iterate    51  f =       1789.5  |proj g|=        1.7958
At iterate    52  f =       1789.4  |proj g|=        1.7966
At iterate    53  f =       1789.4  |proj g|=        1.7892
At iterate    54  f =       1789.2  |proj g|=        1.7826
At iterate    55  f =       1789.2  |proj g|=         1.194
At iterate    56  f =       1789.2  |proj g|=        1.1932
At iterate    57  f =       1789.1  |proj g|=        1.1212
At iterate    58  f =         1789  |proj g|=        1.1135
At iterate    59  f =         1789  |proj g|=         1.875
At iterate    60  f =         1789  |proj g|=        1.7381
At iterate    61  f =       1788.9  |proj g|=        1.4756
At iterate    62  f =       1788.9  |proj g|=        1.7829
At iterate    63  f =       1788.8  |proj g|=        1.7807
At iterate    64  f =       1788.8  |proj g|=        1.7842
At iterate    65  f =       1788.8  |proj g|=       0.55799
At iterate    66  f =       1788.8  |proj g|=       0.54893
At iterate    67  f =       1788.8  |proj g|=         1.162
At iterate    68  f =       1788.8  |proj g|=        1.8723
At iterate    69  f =       1788.8  |proj g|=        1.0113
At iterate    70  f =       1788.8  |proj g|=       0.42081
At iterate    71  f =       1788.8  |proj g|=       0.40831
At iterate    72  f =       1788.8  |proj g|=       0.37755
At iterate    73  f =       1788.7  |proj g|=       0.55045
At iterate    74  f =       1788.7  |proj g|=        1.8706
At iterate    75  f =       1788.7  |proj g|=       0.62034
At iterate    76  f =       1788.7  |proj g|=       0.32525
At iterate    77  f =       1788.7  |proj g|=       0.30756
At iterate    78  f =       1788.7  |proj g|=       0.39267
At iterate    79  f =       1788.7  |proj g|=       0.63073
At iterate    80  f =       1788.7  |proj g|=       0.54849
At iterate    81  f =       1788.7  |proj g|=       0.37813
At iterate    82  f =       1788.7  |proj g|=       0.28969
At iterate    83  f =       1788.7  |proj g|=       0.40478
At iterate    84  f =       1788.7  |proj g|=        0.5533
At iterate    85  f =       1788.7  |proj g|=       0.68968
At iterate    86  f =       1788.7  |proj g|=        1.2444
At iterate    87  f =       1788.7  |proj g|=      0.099835
At iterate    88  f =       1788.7  |proj g|=      0.099618
At iterate    89  f =       1788.7  |proj g|=        0.1305
At iterate    90  f =       1788.7  |proj g|=       0.13061
At iterate    91  f =       1788.7  |proj g|=       0.47715
At iterate    92  f =       1788.7  |proj g|=       0.21897
At iterate    93  f =       1788.7  |proj g|=        0.1916
At iterate    94  f =       1788.7  |proj g|=       0.22821
At iterate    95  f =       1788.7  |proj g|=        1.0767
At iterate    96  f =       1788.7  |proj g|=        0.3527
At iterate    97  f =       1788.7  |proj g|=       0.33189
At iterate    98  f =       1788.7  |proj g|=       0.21486
At iterate    99  f =       1788.7  |proj g|=        0.1876
At iterate   100  f =       1788.7  |proj g|=       0.22777
At iterate   101  f =       1788.7  |proj g|=       0.53975
final  value 1788.683283 
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1844.2  |proj g|=        1.901
At iterate     1  f =       1829.2  |proj g|=        1.7095
At iterate     2  f =       1826.3  |proj g|=        1.7222
At iterate     3  f =       1817.4  |proj g|=        1.8951
At iterate     4  f =       1810.3  |proj g|=        1.9269
At iterate     5  f =       1809.7  |proj g|=        1.3481
At iterate     6  f =       1808.7  |proj g|=        1.3548
At iterate     7  f =       1807.2  |proj g|=        1.9013
At iterate     8  f =       1805.9  |proj g|=        1.9149
At iterate     9  f =       1805.4  |proj g|=        1.8995
At iterate    10  f =       1803.8  |proj g|=        1.1077
At iterate    11  f =       1803.5  |proj g|=        1.9043
At iterate    12  f =       1803.4  |proj g|=        1.0665
At iterate    13  f =       1802.9  |proj g|=        1.9108
At iterate    14  f =       1802.7  |proj g|=        1.0578
At iterate    15  f =       1801.1  |proj g|=        1.9175
At iterate    16  f =       1800.8  |proj g|=        1.6693
At iterate    17  f =       1800.1  |proj g|=        1.6028
At iterate    18  f =       1799.4  |proj g|=        1.9058
At iterate    19  f =         1798  |proj g|=        1.8989
At iterate    20  f =       1797.2  |proj g|=       0.94213
At iterate    21  f =       1796.9  |proj g|=        1.8982
At iterate    22  f =       1796.3  |proj g|=        1.3185
At iterate    23  f =       1794.8  |proj g|=         1.524
At iterate    24  f =       1794.4  |proj g|=        1.5174
At iterate    25  f =       1793.8  |proj g|=        1.8934
At iterate    26  f =       1793.7  |proj g|=        1.8903
At iterate    27  f =       1793.6  |proj g|=        1.4863
At iterate    28  f =       1793.4  |proj g|=        1.4153
At iterate    29  f =       1793.3  |proj g|=        1.6275
At iterate    30  f =       1793.1  |proj g|=        1.7603
At iterate    31  f =       1792.9  |proj g|=        1.8829
At iterate    32  f =       1792.7  |proj g|=        1.8851
At iterate    33  f =       1792.5  |proj g|=        1.7591
At iterate    34  f =       1792.4  |proj g|=       0.74821
At iterate    35  f =         1792  |proj g|=       0.67033
At iterate    36  f =       1791.2  |proj g|=       0.88902
At iterate    37  f =         1791  |proj g|=        1.8859
At iterate    38  f =       1790.9  |proj g|=        0.8918
At iterate    39  f =       1790.9  |proj g|=         1.885
At iterate    40  f =       1790.9  |proj g|=       0.89031
At iterate    41  f =       1790.8  |proj g|=       0.88962
At iterate    42  f =       1790.8  |proj g|=         1.884
At iterate    43  f =       1790.7  |proj g|=        1.6486
At iterate    44  f =       1790.7  |proj g|=       0.87749
At iterate    45  f =       1790.5  |proj g|=       0.87374
At iterate    46  f =       1790.2  |proj g|=       0.92403
At iterate    47  f =         1790  |proj g|=       0.89554
At iterate    48  f =       1789.9  |proj g|=       0.85871
At iterate    49  f =       1789.8  |proj g|=       0.81672
At iterate    50  f =       1789.8  |proj g|=        1.7848
At iterate    51  f =       1789.8  |proj g|=        1.7893
At iterate    52  f =       1789.7  |proj g|=        1.7884
At iterate    53  f =       1789.7  |proj g|=       0.75862
At iterate    54  f =       1789.7  |proj g|=        0.7675
At iterate    55  f =       1789.6  |proj g|=       0.77272
At iterate    56  f =       1789.5  |proj g|=        1.8809
At iterate    57  f =       1789.4  |proj g|=       0.51739
At iterate    58  f =       1789.3  |proj g|=        1.8875
At iterate    59  f =       1789.2  |proj g|=        1.8826
At iterate    60  f =       1789.2  |proj g|=        1.7779
At iterate    61  f =       1789.2  |proj g|=        1.3471
At iterate    62  f =       1789.1  |proj g|=       0.54392
At iterate    63  f =       1789.1  |proj g|=       0.86088
At iterate    64  f =       1789.1  |proj g|=        1.8811
At iterate    65  f =       1789.1  |proj g|=        1.8811
At iterate    66  f =         1789  |proj g|=        1.8798
At iterate    67  f =       1788.9  |proj g|=        1.8743
At iterate    68  f =       1788.9  |proj g|=       0.80123
At iterate    69  f =       1788.9  |proj g|=       0.78927
At iterate    70  f =       1788.8  |proj g|=        1.8759
At iterate    71  f =       1788.8  |proj g|=       0.71153
At iterate    72  f =       1788.7  |proj g|=       0.52376
At iterate    73  f =       1788.7  |proj g|=       0.37546
At iterate    74  f =       1788.7  |proj g|=       0.32805
At iterate    75  f =       1788.7  |proj g|=       0.31573
At iterate    76  f =       1788.7  |proj g|=       0.29298
At iterate    77  f =       1788.7  |proj g|=        1.8724
At iterate    78  f =       1788.7  |proj g|=       0.77072
At iterate    79  f =       1788.7  |proj g|=       0.42465
At iterate    80  f =       1788.7  |proj g|=       0.30258
At iterate    81  f =       1788.7  |proj g|=       0.55782
At iterate    82  f =       1788.7  |proj g|=       0.59316
At iterate    83  f =       1788.7  |proj g|=       0.61489
At iterate    84  f =       1788.7  |proj g|=       0.50715
At iterate    85  f =       1788.7  |proj g|=       0.16702
At iterate    86  f =       1788.7  |proj g|=       0.28302
At iterate    87  f =       1788.7  |proj g|=       0.33463
At iterate    88  f =       1788.7  |proj g|=       0.67044
At iterate    89  f =       1788.7  |proj g|=       0.66924
At iterate    90  f =       1788.7  |proj g|=       0.21373
At iterate    91  f =       1788.7  |proj g|=       0.21343
At iterate    92  f =       1788.7  |proj g|=       0.12906
At iterate    93  f =       1788.7  |proj g|=       0.23661
At iterate    94  f =       1788.7  |proj g|=       0.15096
At iterate    95  f =       1788.7  |proj g|=      0.074849
At iterate    96  f =       1788.7  |proj g|=      0.089774
At iterate    97  f =       1788.7  |proj g|=       0.15525
At iterate    98  f =       1788.7  |proj g|=       0.20395
At iterate    99  f =       1788.7  |proj g|=       0.19495
At iterate   100  f =       1788.7  |proj g|=       0.12942
At iterate   101  f =       1788.7  |proj g|=       0.12935
final  value 1788.682106 
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1844.5  |proj g|=       1.8965
At iterate     1  f =       1837.2  |proj g|=        1.4588
At iterate     2  f =       1830.9  |proj g|=        1.6454
ys=-1.833e-01  -gs= 5.897e+00, BFGS update SKIPPED
At iterate     3  f =       1824.7  |proj g|=        1.8975
At iterate     4  f =       1821.1  |proj g|=        1.8605
ys=-6.464e-01  -gs= 3.354e+00, BFGS update SKIPPED
At iterate     5  f =         1815  |proj g|=        1.9047
At iterate     6  f =       1810.3  |proj g|=        1.8804
At iterate     7  f =       1809.5  |proj g|=        1.5698
At iterate     8  f =       1808.6  |proj g|=        1.8781
At iterate     9  f =         1808  |proj g|=        1.8382
At iterate    10  f =       1807.5  |proj g|=         1.657
At iterate    11  f =       1806.9  |proj g|=        1.8304
At iterate    12  f =       1806.3  |proj g|=        1.8426
At iterate    13  f =       1804.9  |proj g|=        1.8294
At iterate    14  f =       1804.6  |proj g|=        1.6177
At iterate    15  f =       1804.3  |proj g|=        1.5993
At iterate    16  f =       1804.1  |proj g|=         1.853
At iterate    17  f =       1803.7  |proj g|=         1.568
At iterate    18  f =       1803.5  |proj g|=        1.8548
At iterate    19  f =       1798.5  |proj g|=        1.8698
At iterate    20  f =         1798  |proj g|=        1.5165
At iterate    21  f =       1796.5  |proj g|=        1.8676
At iterate    22  f =       1796.4  |proj g|=        1.5568
At iterate    23  f =       1796.2  |proj g|=        1.5642
At iterate    24  f =       1793.6  |proj g|=        1.5362
At iterate    25  f =       1793.3  |proj g|=        1.8045
At iterate    26  f =       1793.1  |proj g|=        1.8804
At iterate    27  f =       1792.9  |proj g|=        1.8841
At iterate    28  f =       1792.6  |proj g|=        1.8834
At iterate    29  f =         1792  |proj g|=        1.4467
At iterate    30  f =       1791.8  |proj g|=        1.2888
At iterate    31  f =       1790.9  |proj g|=        1.8787
At iterate    32  f =       1790.7  |proj g|=        1.0173
At iterate    33  f =       1790.5  |proj g|=       0.98947
At iterate    34  f =       1790.1  |proj g|=        1.8805
At iterate    35  f =         1790  |proj g|=        1.4888
At iterate    36  f =         1790  |proj g|=       0.70684
At iterate    37  f =         1790  |proj g|=       0.91985
At iterate    38  f =         1790  |proj g|=        1.6393
At iterate    39  f =       1789.8  |proj g|=         1.882
At iterate    40  f =       1789.7  |proj g|=         1.881
At iterate    41  f =       1789.6  |proj g|=        1.8776
At iterate    42  f =       1789.6  |proj g|=       0.76916
At iterate    43  f =       1789.5  |proj g|=       0.71971
At iterate    44  f =       1789.5  |proj g|=        1.1675
At iterate    45  f =       1789.5  |proj g|=        1.6044
At iterate    46  f =       1789.4  |proj g|=        1.8782
At iterate    47  f =       1789.2  |proj g|=        1.8816
At iterate    48  f =       1789.2  |proj g|=         1.801
At iterate    49  f =       1789.1  |proj g|=        1.7905
At iterate    50  f =       1789.1  |proj g|=       0.92407
At iterate    51  f =       1789.1  |proj g|=       0.83788
At iterate    52  f =         1789  |proj g|=       0.77994
At iterate    53  f =         1789  |proj g|=       0.60603
At iterate    54  f =       1788.9  |proj g|=        0.5705
At iterate    55  f =       1788.9  |proj g|=        1.8185
At iterate    56  f =       1788.9  |proj g|=        1.0668
At iterate    57  f =       1788.9  |proj g|=        1.0849
At iterate    58  f =       1788.9  |proj g|=       0.99779
At iterate    59  f =       1788.8  |proj g|=       0.30721
At iterate    60  f =       1788.8  |proj g|=       0.24539
At iterate    61  f =       1788.8  |proj g|=        1.0934
At iterate    62  f =       1788.8  |proj g|=        1.1857
At iterate    63  f =       1788.8  |proj g|=       0.40884
At iterate    64  f =       1788.8  |proj g|=       0.35038
At iterate    65  f =       1788.8  |proj g|=       0.57915
At iterate    66  f =       1788.8  |proj g|=       0.55331
At iterate    67  f =       1788.8  |proj g|=       0.75052
At iterate    68  f =       1788.7  |proj g|=       0.99308
At iterate    69  f =       1788.7  |proj g|=        1.7903
At iterate    70  f =       1788.7  |proj g|=       0.57204
At iterate    71  f =       1788.7  |proj g|=       0.30753
At iterate    72  f =       1788.7  |proj g|=       0.40401
At iterate    73  f =       1788.7  |proj g|=        1.4657
At iterate    74  f =       1788.7  |proj g|=       0.55761
At iterate    75  f =       1788.7  |proj g|=       0.49068
At iterate    76  f =       1788.7  |proj g|=       0.42876
At iterate    77  f =       1788.7  |proj g|=       0.46333
At iterate    78  f =       1788.7  |proj g|=         1.208
At iterate    79  f =       1788.7  |proj g|=        0.6762
At iterate    80  f =       1788.7  |proj g|=       0.49907
At iterate    81  f =       1788.7  |proj g|=       0.33404
At iterate    82  f =       1788.7  |proj g|=       0.16195
At iterate    83  f =       1788.7  |proj g|=       0.40693
At iterate    84  f =       1788.7  |proj g|=       0.19521
At iterate    85  f =       1788.7  |proj g|=       0.16413
At iterate    86  f =       1788.7  |proj g|=       0.13068
At iterate    87  f =       1788.7  |proj g|=      0.061973
At iterate    88  f =       1788.7  |proj g|=        0.1542
At iterate    89  f =       1788.7  |proj g|=       0.18834
At iterate    90  f =       1788.7  |proj g|=       0.11761
At iterate    91  f =       1788.7  |proj g|=       0.11372
At iterate    92  f =       1788.7  |proj g|=       0.13048
At iterate    93  f =       1788.7  |proj g|=       0.41079
At iterate    94  f =       1788.7  |proj g|=       0.13047
At iterate    95  f =       1788.7  |proj g|=      0.060687
At iterate    96  f =       1788.7  |proj g|=      0.061571
At iterate    97  f =       1788.7  |proj g|=       0.13531
At iterate    98  f =       1788.7  |proj g|=       0.20686
At iterate    99  f =       1788.7  |proj g|=       0.13803
At iterate   100  f =       1788.7  |proj g|=       0.44702
At iterate   101  f =       1788.7  |proj g|=       0.17965
final  value 1788.680770 
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1844.5  |proj g|=       1.9627
At iterate     1  f =       1841.8  |proj g|=        1.5053
At iterate     2  f =       1835.2  |proj g|=        1.5189
At iterate     3  f =       1832.8  |proj g|=        1.9281
At iterate     4  f =         1828  |proj g|=        1.2229
At iterate     5  f =       1823.6  |proj g|=        1.0775
At iterate     6  f =         1822  |proj g|=        1.9408
At iterate     7  f =       1809.7  |proj g|=        1.8918
At iterate     8  f =         1803  |proj g|=         1.845
At iterate     9  f =       1802.6  |proj g|=        1.7398
At iterate    10  f =       1801.4  |proj g|=        1.2403
At iterate    11  f =       1799.9  |proj g|=        1.8317
At iterate    12  f =         1799  |proj g|=        1.8231
At iterate    13  f =       1797.6  |proj g|=       0.99073
At iterate    14  f =       1795.2  |proj g|=        1.0283
At iterate    15  f =       1793.4  |proj g|=        1.8782
At iterate    16  f =       1792.8  |proj g|=        1.7883
At iterate    17  f =       1792.7  |proj g|=         1.544
At iterate    18  f =       1792.6  |proj g|=        1.8611
At iterate    19  f =       1792.5  |proj g|=        1.8626
At iterate    20  f =       1791.9  |proj g|=         1.867
At iterate    21  f =       1791.7  |proj g|=        1.4722
At iterate    22  f =       1791.6  |proj g|=         1.456
At iterate    23  f =       1791.3  |proj g|=        1.7959
At iterate    24  f =         1791  |proj g|=        1.7963
At iterate    25  f =       1790.9  |proj g|=        1.8665
At iterate    26  f =       1790.7  |proj g|=        1.8672
At iterate    27  f =       1790.7  |proj g|=       0.78433
At iterate    28  f =       1790.6  |proj g|=       0.78289
At iterate    29  f =       1790.6  |proj g|=        0.7744
At iterate    30  f =       1790.4  |proj g|=        1.8684
At iterate    31  f =       1790.4  |proj g|=        1.7928
At iterate    32  f =       1790.3  |proj g|=       0.80236
At iterate    33  f =       1790.3  |proj g|=       0.82648
At iterate    34  f =       1790.1  |proj g|=       0.87867
At iterate    35  f =         1790  |proj g|=        0.8733
At iterate    36  f =       1789.8  |proj g|=        1.7905
At iterate    37  f =       1789.7  |proj g|=        1.8743
At iterate    38  f =       1789.7  |proj g|=        1.8715
At iterate    39  f =       1789.6  |proj g|=       0.80037
At iterate    40  f =       1789.6  |proj g|=       0.80415
At iterate    41  f =       1789.5  |proj g|=       0.80189
At iterate    42  f =       1789.5  |proj g|=       0.78967
At iterate    43  f =       1789.4  |proj g|=       0.75391
At iterate    44  f =       1789.3  |proj g|=         1.872
At iterate    45  f =       1789.3  |proj g|=       0.68717
At iterate    46  f =       1789.1  |proj g|=        1.4814
At iterate    47  f =       1789.1  |proj g|=         1.487
At iterate    48  f =         1789  |proj g|=        1.8742
At iterate    49  f =         1789  |proj g|=        1.4876
At iterate    50  f =       1788.9  |proj g|=        1.7243
At iterate    51  f =       1788.9  |proj g|=       0.49094
At iterate    52  f =       1788.9  |proj g|=        0.5966
At iterate    53  f =       1788.9  |proj g|=       0.59535
At iterate    54  f =       1788.9  |proj g|=        1.8733
At iterate    55  f =       1788.8  |proj g|=       0.93604
At iterate    56  f =       1788.8  |proj g|=       0.85211
At iterate    57  f =       1788.8  |proj g|=        1.4377
At iterate    58  f =       1788.8  |proj g|=        1.4772
At iterate    59  f =       1788.8  |proj g|=       0.78361
At iterate    60  f =       1788.8  |proj g|=        1.8383
At iterate    61  f =       1788.8  |proj g|=        0.3324
At iterate    62  f =       1788.8  |proj g|=       0.29873
At iterate    63  f =       1788.8  |proj g|=       0.94247
At iterate    64  f =       1788.8  |proj g|=       0.30321
At iterate    65  f =       1788.8  |proj g|=       0.32765
At iterate    66  f =       1788.7  |proj g|=       0.48866
At iterate    67  f =       1788.7  |proj g|=        1.1055
At iterate    68  f =       1788.7  |proj g|=       0.71111
At iterate    69  f =       1788.7  |proj g|=       0.50294
At iterate    70  f =       1788.7  |proj g|=       0.95867
At iterate    71  f =       1788.7  |proj g|=       0.49225
At iterate    72  f =       1788.7  |proj g|=       0.33465
At iterate    73  f =       1788.7  |proj g|=        0.3246
At iterate    74  f =       1788.7  |proj g|=       0.27462
At iterate    75  f =       1788.7  |proj g|=       0.22952
At iterate    76  f =       1788.7  |proj g|=       0.47635
At iterate    77  f =       1788.7  |proj g|=       0.36506
At iterate    78  f =       1788.7  |proj g|=        0.2705
At iterate    79  f =       1788.7  |proj g|=       0.35958
At iterate    80  f =       1788.7  |proj g|=       0.21413
At iterate    81  f =       1788.7  |proj g|=      0.083278
At iterate    82  f =       1788.7  |proj g|=      0.069125
At iterate    83  f =       1788.7  |proj g|=      0.044831
At iterate    84  f =       1788.7  |proj g|=       0.13009
At iterate    85  f =       1788.7  |proj g|=      0.083972
At iterate    86  f =       1788.7  |proj g|=      0.038698
At iterate    87  f =       1788.7  |proj g|=      0.055152
At iterate    88  f =       1788.7  |proj g|=       0.02741
At iterate    89  f =       1788.7  |proj g|=      0.026779
At iterate    90  f =       1788.7  |proj g|=        0.0423
At iterate    91  f =       1788.7  |proj g|=      0.065137
At iterate    92  f =       1788.7  |proj g|=      0.046154
At iterate    93  f =       1788.7  |proj g|=      0.041203
At iterate    94  f =       1788.7  |proj g|=      0.058948
At iterate    95  f =       1788.7  |proj g|=      0.015991
At iterate    96  f =       1788.7  |proj g|=      0.021177

iterations 96
function evaluations 123
segments explored during Cauchy searches 101
BFGS updates skipped 0
active bounds at final generalized Cauchy point 23
norm of the final projected gradient 0.0211773
final function value 1788.68

F = 1788.68
final  value 1788.679245 
converged

* The 4 best values (multistart) obtained are:
 1788.683 1788.682 1788.681 1788.679 
* The model corresponding to the best one (1788.679) is stored. 

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2 
  - best initial criterion value(s) :  -6193.86 -6193.966 -6195.06 -6196.14 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       6193.9  |proj g|=       1.9414
At iterate     1  f =       6099.8  |proj g|=        1.5695
At iterate     2  f =       6091.5  |proj g|=         1.903
At iterate     3  f =       6089.3  |proj g|=        1.8889
At iterate     4  f =       6083.7  |proj g|=        1.8194
At iterate     5  f =       6082.8  |proj g|=        1.8071
At iterate     6  f =       6079.7  |proj g|=        1.9016
At iterate     7  f =       6077.9  |proj g|=        1.8974
At iterate     8  f =       6068.1  |proj g|=        1.4921
At iterate     9  f =       6066.4  |proj g|=        1.4399
At iterate    10  f =       6064.5  |proj g|=        1.9119
At iterate    11  f =       6063.5  |proj g|=        1.5745
At iterate    12  f =       6060.4  |proj g|=        1.5219
At iterate    13  f =         6060  |proj g|=        1.5145
At iterate    14  f =         6059  |proj g|=        1.9074
At iterate    15  f =       6057.9  |proj g|=        1.8963
At iterate    16  f =       6056.2  |proj g|=        1.2182
At iterate    17  f =       6054.7  |proj g|=        1.8811
At iterate    18  f =       6053.9  |proj g|=        1.1028
At iterate    19  f =       6053.2  |proj g|=        1.0823
At iterate    20  f =       6053.1  |proj g|=         1.888
At iterate    21  f =       6052.9  |proj g|=        1.0711
At iterate    22  f =       6052.7  |proj g|=        1.4231
At iterate    23  f =       6052.5  |proj g|=        1.8888
At iterate    24  f =       6052.4  |proj g|=        1.3347
At iterate    25  f =       6052.4  |proj g|=        1.0392
At iterate    26  f =       6052.4  |proj g|=        1.8869
At iterate    27  f =       6052.3  |proj g|=        1.0399
At iterate    28  f =       6052.3  |proj g|=        1.0381
At iterate    29  f =       6052.3  |proj g|=         1.023
At iterate    30  f =       6052.1  |proj g|=       0.80119
At iterate    31  f =       6051.9  |proj g|=       0.74419
At iterate    32  f =       6051.9  |proj g|=       0.69194
At iterate    33  f =       6051.8  |proj g|=        1.7041
At iterate    34  f =       6051.8  |proj g|=       0.83632
At iterate    35  f =       6051.8  |proj g|=       0.83293
At iterate    36  f =       6051.7  |proj g|=       0.64953
At iterate    37  f =       6051.7  |proj g|=       0.48825
At iterate    38  f =       6051.7  |proj g|=       0.28695
At iterate    39  f =       6051.7  |proj g|=        1.6661
At iterate    40  f =       6051.7  |proj g|=       0.46224
At iterate    41  f =       6051.7  |proj g|=       0.15696
At iterate    42  f =       6051.7  |proj g|=       0.14672
At iterate    43  f =       6051.7  |proj g|=       0.11677
At iterate    44  f =       6051.7  |proj g|=      0.048306
At iterate    45  f =       6051.7  |proj g|=      0.048253
At iterate    46  f =       6051.7  |proj g|=       0.09721
At iterate    47  f =       6051.7  |proj g|=       0.11672
At iterate    48  f =       6051.7  |proj g|=         0.103
At iterate    49  f =       6051.7  |proj g|=       0.10676
At iterate    50  f =       6051.7  |proj g|=      0.012401

iterations 50
function evaluations 56
segments explored during Cauchy searches 62
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0124006
final function value 6051.72

F = 6051.72
final  value 6051.720344 
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=         6194  |proj g|=        1.896
At iterate     1  f =       6096.2  |proj g|=        1.9139
At iterate     2  f =       6095.7  |proj g|=        1.9117
At iterate     3  f =       6092.7  |proj g|=        1.9101
At iterate     4  f =       6092.3  |proj g|=        1.9112
At iterate     5  f =       6091.6  |proj g|=        1.9106
ys=-1.660e-03  -gs= 7.046e-01, BFGS update SKIPPED
At iterate     6  f =       6088.6  |proj g|=        1.9268
At iterate     7  f =         6080  |proj g|=        1.9345
ys=-1.615e-02  -gs= 8.492e+00, BFGS update SKIPPED
At iterate     8  f =         6077  |proj g|=        1.9229
At iterate     9  f =       6069.1  |proj g|=        1.8154
At iterate    10  f =       6068.5  |proj g|=        1.6384
At iterate    11  f =         6068  |proj g|=        1.9184
At iterate    12  f =       6067.9  |proj g|=        1.9211
At iterate    13  f =         6067  |proj g|=        1.9248
At iterate    14  f =       6063.9  |proj g|=        1.9201
At iterate    15  f =       6062.4  |proj g|=        1.9124
At iterate    16  f =       6061.1  |proj g|=        1.7501
At iterate    17  f =       6059.5  |proj g|=        1.9052
At iterate    18  f =       6059.1  |proj g|=        1.7219
At iterate    19  f =       6057.7  |proj g|=        1.7267
At iterate    20  f =       6056.9  |proj g|=        1.7136
At iterate    21  f =       6056.1  |proj g|=        1.6472
At iterate    22  f =       6055.3  |proj g|=        1.8976
At iterate    23  f =       6055.2  |proj g|=        1.6393
At iterate    24  f =         6055  |proj g|=         1.629
At iterate    25  f =       6054.6  |proj g|=        1.8972
At iterate    26  f =       6054.6  |proj g|=        1.5981
At iterate    27  f =       6054.5  |proj g|=        1.5928
At iterate    28  f =       6054.4  |proj g|=        1.5774
At iterate    29  f =       6054.2  |proj g|=        1.8957
At iterate    30  f =       6054.1  |proj g|=        1.8949
At iterate    31  f =       6052.6  |proj g|=        1.2987
At iterate    32  f =       6052.3  |proj g|=        1.2399
At iterate    33  f =       6051.9  |proj g|=       0.96675
At iterate    34  f =       6051.8  |proj g|=       0.96791
At iterate    35  f =       6051.8  |proj g|=       0.94212
At iterate    36  f =       6051.8  |proj g|=        1.8852
At iterate    37  f =       6051.8  |proj g|=       0.71337
At iterate    38  f =       6051.8  |proj g|=       0.71134
At iterate    39  f =       6051.7  |proj g|=        0.3276
At iterate    40  f =       6051.7  |proj g|=       0.19604
At iterate    41  f =       6051.7  |proj g|=       0.14445
At iterate    42  f =       6051.7  |proj g|=       0.13129
At iterate    43  f =       6051.7  |proj g|=       0.11073
At iterate    44  f =       6051.7  |proj g|=       0.65001
At iterate    45  f =       6051.7  |proj g|=       0.50626
At iterate    46  f =       6051.7  |proj g|=       0.21399
At iterate    47  f =       6051.7  |proj g|=        0.1184
At iterate    48  f =       6051.7  |proj g|=       0.11665
At iterate    49  f =       6051.7  |proj g|=      0.076832
At iterate    50  f =       6051.7  |proj g|=       0.11684
At iterate    51  f =       6051.7  |proj g|=        0.1167
At iterate    52  f =       6051.7  |proj g|=       0.09314
At iterate    53  f =       6051.7  |proj g|=       0.10182
At iterate    54  f =       6051.7  |proj g|=      0.050777
At iterate    55  f =       6051.7  |proj g|=      0.025595
At iterate    56  f =       6051.7  |proj g|=       0.11416

iterations 56
function evaluations 73
segments explored during Cauchy searches 68
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.114159
final function value 6051.72

F = 6051.72
final  value 6051.720361 
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       6195.1  |proj g|=       1.8444
At iterate     1  f =       6190.4  |proj g|=         1.567
At iterate     2  f =       6180.2  |proj g|=        1.8736
At iterate     3  f =       6172.8  |proj g|=        1.8258
ys=-8.212e-01  -gs= 6.952e+00, BFGS update SKIPPED
At iterate     4  f =       6153.1  |proj g|=        1.7389
ys=-2.505e+00  -gs= 1.744e+01, BFGS update SKIPPED
At iterate     5  f =       6141.1  |proj g|=         1.694
At iterate     6  f =       6130.4  |proj g|=        1.9376
At iterate     7  f =         6113  |proj g|=        1.8015
At iterate     8  f =       6107.7  |proj g|=        1.8026
At iterate     9  f =       6085.2  |proj g|=        1.6781
At iterate    10  f =         6084  |proj g|=        1.8835
At iterate    11  f =       6082.6  |proj g|=        1.8806
At iterate    12  f =       6072.2  |proj g|=        1.8738
At iterate    13  f =       6070.1  |proj g|=        1.4886
At iterate    14  f =       6065.8  |proj g|=         1.381
At iterate    15  f =       6064.4  |proj g|=        1.3835
At iterate    16  f =         6064  |proj g|=        1.8899
At iterate    17  f =       6062.2  |proj g|=        1.8814
At iterate    18  f =       6061.5  |proj g|=         1.421
At iterate    19  f =       6059.1  |proj g|=        1.4174
At iterate    20  f =       6058.3  |proj g|=        1.4164
At iterate    21  f =       6056.7  |proj g|=        1.4007
At iterate    22  f =       6055.7  |proj g|=         1.901
At iterate    23  f =       6055.3  |proj g|=        1.4749
At iterate    24  f =       6054.9  |proj g|=        1.4674
At iterate    25  f =       6054.7  |proj g|=        1.8948
At iterate    26  f =       6054.5  |proj g|=        1.8955
At iterate    27  f =       6054.4  |proj g|=        1.8936
At iterate    28  f =       6053.8  |proj g|=        1.1737
At iterate    29  f =       6053.3  |proj g|=        1.8945
At iterate    30  f =       6053.2  |proj g|=        1.8917
At iterate    31  f =       6053.1  |proj g|=         1.353
At iterate    32  f =       6052.9  |proj g|=        1.4434
At iterate    33  f =       6052.8  |proj g|=        1.8925
At iterate    34  f =       6052.7  |proj g|=         1.441
At iterate    35  f =       6052.4  |proj g|=        1.4288
At iterate    36  f =       6052.3  |proj g|=        1.0089
At iterate    37  f =       6052.2  |proj g|=        1.8889
At iterate    38  f =       6052.1  |proj g|=        1.8888
At iterate    39  f =       6052.1  |proj g|=        1.6079
At iterate    40  f =         6052  |proj g|=        1.2944
At iterate    41  f =         6052  |proj g|=        1.1205
At iterate    42  f =         6052  |proj g|=        1.8863
At iterate    43  f =         6052  |proj g|=         1.429
At iterate    44  f =       6051.9  |proj g|=       0.92281
At iterate    45  f =       6051.9  |proj g|=        1.3317
At iterate    46  f =       6051.9  |proj g|=        1.3044
At iterate    47  f =       6051.9  |proj g|=       0.78826
At iterate    48  f =       6051.9  |proj g|=       0.80595
At iterate    49  f =       6051.9  |proj g|=       0.91726
At iterate    50  f =       6051.8  |proj g|=        1.8835
At iterate    51  f =       6051.7  |proj g|=        1.8839
At iterate    52  f =       6051.7  |proj g|=       0.62815
At iterate    53  f =       6051.7  |proj g|=       0.59062
At iterate    54  f =       6051.7  |proj g|=       0.09207
At iterate    55  f =       6051.7  |proj g|=      0.069011
At iterate    56  f =       6051.7  |proj g|=      0.016776
At iterate    57  f =       6051.7  |proj g|=      0.007937

iterations 57
function evaluations 64
segments explored during Cauchy searches 65
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.00793695
final function value 6051.72

F = 6051.72
final  value 6051.720335 
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       6196.1  |proj g|=       1.9441
At iterate     1  f =       6158.6  |proj g|=        1.8479
At iterate     2  f =       6131.1  |proj g|=        1.8288
ys=-4.125e+00  -gs= 2.537e+01, BFGS update SKIPPED
At iterate     3  f =       6093.3  |proj g|=        1.7964
At iterate     4  f =       6091.6  |proj g|=        1.8966
At iterate     5  f =       6089.7  |proj g|=        1.7487
At iterate     6  f =       6078.4  |proj g|=        1.7556
At iterate     7  f =       6069.5  |proj g|=        1.7184
At iterate     8  f =       6066.2  |proj g|=        1.6136
At iterate     9  f =       6065.3  |proj g|=        1.9357
At iterate    10  f =       6063.3  |proj g|=        1.9207
At iterate    11  f =       6061.3  |proj g|=         1.909
At iterate    12  f =       6059.9  |proj g|=        1.8995
At iterate    13  f =       6057.7  |proj g|=        1.7114
At iterate    14  f =       6057.6  |proj g|=        1.7106
At iterate    15  f =       6057.1  |proj g|=        1.9015
At iterate    16  f =       6056.4  |proj g|=        1.9085
At iterate    17  f =         6056  |proj g|=        1.8998
At iterate    18  f =       6055.9  |proj g|=        1.6836
At iterate    19  f =       6055.7  |proj g|=        1.6795
At iterate    20  f =       6055.6  |proj g|=        1.8978
At iterate    21  f =       6055.5  |proj g|=        1.8959
At iterate    22  f =       6055.3  |proj g|=        1.6477
At iterate    23  f =         6055  |proj g|=        1.6138
At iterate    24  f =       6054.5  |proj g|=        1.5619
At iterate    25  f =       6054.1  |proj g|=        1.8977
At iterate    26  f =         6054  |proj g|=        1.8975
At iterate    27  f =       6053.4  |proj g|=        1.4395
At iterate    28  f =         6053  |proj g|=         1.889
At iterate    29  f =       6052.3  |proj g|=        1.8807
At iterate    30  f =         6052  |proj g|=       0.72993
At iterate    31  f =       6051.9  |proj g|=       0.84726
At iterate    32  f =       6051.8  |proj g|=        1.8843
At iterate    33  f =       6051.8  |proj g|=       0.47819
At iterate    34  f =       6051.8  |proj g|=        1.8846
At iterate    35  f =       6051.7  |proj g|=       0.70317
At iterate    36  f =       6051.7  |proj g|=       0.58719
At iterate    37  f =       6051.7  |proj g|=       0.53426
At iterate    38  f =       6051.7  |proj g|=       0.28181
At iterate    39  f =       6051.7  |proj g|=       0.27471
At iterate    40  f =       6051.7  |proj g|=       0.71023
At iterate    41  f =       6051.7  |proj g|=        0.7442
At iterate    42  f =       6051.7  |proj g|=       0.17579
At iterate    43  f =       6051.7  |proj g|=       0.12481
At iterate    44  f =       6051.7  |proj g|=       0.11666
At iterate    45  f =       6051.7  |proj g|=       0.80637
At iterate    46  f =       6051.7  |proj g|=      0.047588
At iterate    47  f =       6051.7  |proj g|=      0.041069
At iterate    48  f =       6051.7  |proj g|=        0.1168
At iterate    49  f =       6051.7  |proj g|=       0.11675
At iterate    50  f =       6051.7  |proj g|=      0.043049
At iterate    51  f =       6051.7  |proj g|=      0.042103
At iterate    52  f =       6051.7  |proj g|=      0.065384
At iterate    53  f =       6051.7  |proj g|=       0.11522
At iterate    54  f =       6051.7  |proj g|=      0.044015

iterations 54
function evaluations 68
segments explored during Cauchy searches 64
BFGS updates skipped 1
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0440146
final function value 6051.72

F = 6051.72
final  value 6051.720388 
converged

* The 4 best values (multistart) obtained are:
 6051.72 6051.72 6051.72 6051.72 
* The model corresponding to the best one (6051.72) is stored. 

[1] 10.81038
[1] 369550.6
[1] 8.353776
[1] 5.220835

1.15 Constrain first and then do sensitivity analysis?

The problem with doing constraint first is that you end up with a non-hypercube shaped input space (corners have been knocked off by constraint), which is not ideal. We might therefore want different sensitivity measures for pre- and post-constrained ensemble.

[1] 8.353776
[1] 8.353615
[1] 5.220835
[1] 5.220903

1.16 Initial sensitivity Matrix

It looks as though bwl_io is very influential across a number of variables, even though it didn’t appear that interesting in the parginal plots. Why is that?

1.17 Re-examine bwl_io

A closer look at bwl_io now that the impact of f0_io has been removed shows a large number of “zero” NPP when bwl_io is at low values, which could well be the source of apparent sensitivity.

Take only higher values of bwl_io for the next round of constraints.

2 Remove ensemble members with low b_wl and see if the sensitivity analyses change and the emulators get better.

# Plot the regular km emulator.

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -5137.861 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       5137.9  |proj g|=       1.8356
At iterate     1  f =       5135.4  |proj g|=        1.7202
At iterate     2  f =       5131.9  |proj g|=           1.6
At iterate     3  f =       5131.2  |proj g|=        1.6826
At iterate     4  f =       5122.2  |proj g|=        1.6838
At iterate     5  f =       5119.7  |proj g|=        1.6697
At iterate     6  f =       5110.5  |proj g|=        1.6667
At iterate     7  f =       5105.6  |proj g|=        1.3939
At iterate     8  f =       5098.3  |proj g|=        1.4145
At iterate     9  f =       5093.3  |proj g|=         1.592
At iterate    10  f =       5091.8  |proj g|=        1.5626
At iterate    11  f =       5089.7  |proj g|=        1.4986
At iterate    12  f =       5088.5  |proj g|=        1.3214
At iterate    13  f =       5085.1  |proj g|=         1.501
At iterate    14  f =       5083.9  |proj g|=        1.4621
At iterate    15  f =       5082.1  |proj g|=        1.4517
At iterate    16  f =       5081.1  |proj g|=         1.494
At iterate    17  f =       5080.6  |proj g|=       0.95799
At iterate    18  f =       5080.4  |proj g|=        1.3623
At iterate    19  f =         5080  |proj g|=       0.56793
At iterate    20  f =       5079.9  |proj g|=        1.4786
At iterate    21  f =       5079.8  |proj g|=       0.46896
At iterate    22  f =       5079.8  |proj g|=       0.28689
At iterate    23  f =       5079.8  |proj g|=       0.57427
At iterate    24  f =       5079.7  |proj g|=       0.57268
At iterate    25  f =       5079.7  |proj g|=       0.36422
At iterate    26  f =       5079.7  |proj g|=       0.16654
At iterate    27  f =       5079.7  |proj g|=       0.22658
At iterate    28  f =       5079.7  |proj g|=       0.27311
At iterate    29  f =       5079.7  |proj g|=       0.21159
At iterate    30  f =       5079.7  |proj g|=       0.52762
At iterate    31  f =       5079.7  |proj g|=        0.2638
At iterate    32  f =       5079.7  |proj g|=       0.15332
At iterate    33  f =       5079.7  |proj g|=       0.16807
At iterate    34  f =       5079.7  |proj g|=        0.2156
At iterate    35  f =       5079.7  |proj g|=      0.093168
At iterate    36  f =       5079.7  |proj g|=      0.069853
At iterate    37  f =       5079.7  |proj g|=       0.11439
At iterate    38  f =       5079.7  |proj g|=       0.18854
At iterate    39  f =       5079.7  |proj g|=       0.14403
At iterate    40  f =       5079.6  |proj g|=      0.075684
At iterate    41  f =       5079.6  |proj g|=      0.045529
At iterate    42  f =       5079.6  |proj g|=      0.029163
At iterate    43  f =       5079.6  |proj g|=      0.028335
At iterate    44  f =       5079.6  |proj g|=      0.081292
At iterate    45  f =       5079.6  |proj g|=       0.03851
At iterate    46  f =       5079.6  |proj g|=      0.011842
At iterate    47  f =       5079.6  |proj g|=      0.010565

iterations 47
function evaluations 54
segments explored during Cauchy searches 60
BFGS updates skipped 0
active bounds at final generalized Cauchy point 20
norm of the final projected gradient 0.0105653
final function value 5079.65

F = 5079.65
final  value 5079.648744 
converged

2.1 With b_wl truncated, we can clearly see the sensitivities change.

2.2 Plot both sensitivity matrices on top of one another.

It looks as though both the absolute value and the change over time are controlled by the same parameters.

2.3 How good are the emulators for each output?

The emulators appear to be at least capturing the broad response for all of the output variables.

First, plot the straight kriging emulators

3 km emulators for change in variables over time

[1] 5.433077
[1] 5.142409
[1] 5.2209
[1] 4.193005
[1] 4.219396
[1] 4.531841
[1] 5.629425
[1] 8.565018
[1] 5.224583
[1] 9.223968
[1] 7.089
[1] 7.78196
[1] 6.900971
[1] 6.929568

Next, plot the twostep glmnet/km emulators

3.1 TwoStep emulator performance for modern values

4 Compare Leave-one-out statistics for the two types of emulator

We use the leave-one-out Mean Absolute Error, expressed as a percentage of the range of the output across the ensemble. We find that the twostep emulatorisn’t significantly more accurate, and is indeed less accurate for tree fraction.

[1] 5.410355
[1] 4.733929
[1] 5.215455
[1] 4.173501
[1] 4.03861
[1] 4.597044
[1] 5.501413
[1] 8.563912
[1] 5.213222
[1] 11.36184
[1] 7.260354
[1] 7.707121
[1] 6.725096
[1] 6.865819

5 History matching and generation of new ensemble members

To use History Matching, we need to specify targets for various model outputs. We treat these targets as “observations” with an uncertainty where a model run marked as “implausible” (beyond 3sd) matches the hard boundaries previously identified by A. Wiltshire as being desirable.

Choose the centre of the (implied) uniform distribution. cs_gb.target = (3000 - 750) / 2 = 1125 cv.target = (800 - 300) / 2 = 250 npp_n_gb.target = (80 - 35) / 2 = 22.5

nbp.target = 0 (gpp.target = 75) (runoff.target = 1)

(to do: visualise implausibility of design and loo emulated values of design as histogram)

5.1 visualise the constraints

5.2 Timeseries that match the constraints.

Plot the original ensemble, plus the stuff that matches the constraints. Include anomalies.

We find that adding a constraint to level 2 constrains the behaviour of the soil carbon pool, but not the vegetation carbon pool. Probably not a surprise, given the Soil Carbon pool (and changes) are much larger than the vegetation carbon pool.

5.3 Total land carbon change

We plot the JULES total land carbon change from 1850, together with “observations” from the Global Carbon Budget (Friedligstein et al, 2020). The “observations are the cumulative sum of the ‘land sink’ column on the”historical observations tab of the data from the paper. It states “The land sink is the average of several dynamic global vegetation models that reproduce the observed mean total land sink of the 1990s.”

5.3.1 Land carbon Sink?

I think to make a fair comparison, we might need to take the diff of the anomaly (or just the timeseries), as the above represents the total change since 1850, and I think all the “land sink” stuff looks at the year-to-year changes.

5.3.2 Identify parts of parameter space where the land carbon sink is at its maximum

What do these parts of parameter space tell us about what we need to do to have the model reproduce reality?

5.3.3 what parts of parameter space are the closest?

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -106.4158 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       106.42  |proj g|=       1.8453
At iterate     1  f =       103.09  |proj g|=         1.819
At iterate     2  f =       97.125  |proj g|=        1.6466
At iterate     3  f =       95.132  |proj g|=        1.6743
At iterate     4  f =       91.572  |proj g|=        1.6622
At iterate     5  f =       89.835  |proj g|=        1.6389
At iterate     6  f =       88.807  |proj g|=        1.5108
At iterate     7  f =       87.513  |proj g|=        1.7187
At iterate     8  f =       86.326  |proj g|=         1.638
At iterate     9  f =       85.522  |proj g|=        1.5586
At iterate    10  f =        85.21  |proj g|=         1.342
At iterate    11  f =        84.63  |proj g|=        1.2524
At iterate    12  f =       84.111  |proj g|=        1.2971
At iterate    13  f =        83.72  |proj g|=        1.1355
At iterate    14  f =        83.52  |proj g|=        1.0089
At iterate    15  f =       82.998  |proj g|=       0.69506
At iterate    16  f =       82.939  |proj g|=        1.7322
At iterate    17  f =       82.779  |proj g|=        1.7049
At iterate    18  f =       82.708  |proj g|=       0.80111
At iterate    19  f =       82.614  |proj g|=       0.71523
At iterate    20  f =       82.575  |proj g|=       0.65009
At iterate    21  f =       82.517  |proj g|=       0.61643
At iterate    22  f =       82.482  |proj g|=       0.53998
At iterate    23  f =        82.35  |proj g|=        1.4933
At iterate    24  f =        82.31  |proj g|=        1.4399
At iterate    25  f =       82.247  |proj g|=       0.83833
At iterate    26  f =       82.179  |proj g|=        0.6743
At iterate    27  f =       82.125  |proj g|=       0.61099
At iterate    28  f =       82.077  |proj g|=       0.69756
At iterate    29  f =       82.068  |proj g|=       0.50643
At iterate    30  f =       82.045  |proj g|=       0.53157
At iterate    31  f =       81.976  |proj g|=        1.2632
At iterate    32  f =       81.806  |proj g|=        1.4306
At iterate    33  f =       81.751  |proj g|=        1.1701
At iterate    34  f =       81.628  |proj g|=        1.1438
At iterate    35  f =       81.556  |proj g|=        1.6474
At iterate    36  f =       81.483  |proj g|=        1.4191
At iterate    37  f =       81.399  |proj g|=       0.67894
At iterate    38  f =       81.324  |proj g|=       0.33289
At iterate    39  f =       81.314  |proj g|=        1.0959
At iterate    40  f =       81.302  |proj g|=       0.32766
At iterate    41  f =       81.295  |proj g|=        0.3297
At iterate    42  f =       81.277  |proj g|=       0.36728
At iterate    43  f =       81.265  |proj g|=       0.30341
At iterate    44  f =       81.247  |proj g|=       0.69733
At iterate    45  f =       81.216  |proj g|=        1.1062
At iterate    46  f =       81.177  |proj g|=        1.1638
At iterate    47  f =       81.155  |proj g|=       0.71712
At iterate    48  f =       81.109  |proj g|=       0.51864
At iterate    49  f =       81.077  |proj g|=       0.41088
At iterate    50  f =       81.059  |proj g|=       0.39031
At iterate    51  f =       81.055  |proj g|=       0.33235
At iterate    52  f =       81.052  |proj g|=       0.14717
At iterate    53  f =       81.049  |proj g|=        0.1497
At iterate    54  f =       81.043  |proj g|=       0.15788
At iterate    55  f =       81.035  |proj g|=       0.29801
At iterate    56  f =       81.033  |proj g|=       0.57956
At iterate    57  f =       81.022  |proj g|=       0.39169
At iterate    58  f =       81.018  |proj g|=       0.16287
At iterate    59  f =       81.015  |proj g|=       0.16427
At iterate    60  f =       81.003  |proj g|=       0.17202
At iterate    61  f =       80.998  |proj g|=       0.55788
At iterate    62  f =       80.983  |proj g|=        0.3853
At iterate    63  f =       80.978  |proj g|=       0.14235
At iterate    64  f =       80.973  |proj g|=       0.16257
At iterate    65  f =       80.969  |proj g|=       0.55124
At iterate    66  f =       80.966  |proj g|=       0.17098
At iterate    67  f =       80.964  |proj g|=       0.15834
At iterate    68  f =       80.961  |proj g|=       0.19492
At iterate    69  f =        80.96  |proj g|=      0.085999
At iterate    70  f =       80.959  |proj g|=       0.20316
At iterate    71  f =       80.957  |proj g|=       0.22804
At iterate    72  f =       80.956  |proj g|=       0.18647
At iterate    73  f =       80.955  |proj g|=      0.058348
At iterate    74  f =       80.955  |proj g|=      0.054364
At iterate    75  f =       80.953  |proj g|=      0.091177
At iterate    76  f =       80.953  |proj g|=       0.31904
At iterate    77  f =       80.951  |proj g|=       0.21525
At iterate    78  f =       80.947  |proj g|=       0.12164
At iterate    79  f =       80.943  |proj g|=       0.12463
At iterate    80  f =       80.941  |proj g|=       0.18732
At iterate    81  f =       80.939  |proj g|=       0.08267
At iterate    82  f =       80.938  |proj g|=       0.24049
At iterate    83  f =       80.938  |proj g|=       0.17522
At iterate    84  f =       80.937  |proj g|=       0.10918
At iterate    85  f =       80.937  |proj g|=      0.079848
At iterate    86  f =       80.936  |proj g|=      0.045969
At iterate    87  f =       80.936  |proj g|=      0.043919
At iterate    88  f =       80.936  |proj g|=       0.12858
At iterate    89  f =       80.936  |proj g|=       0.04115
At iterate    90  f =       80.936  |proj g|=      0.025463
At iterate    91  f =       80.936  |proj g|=      0.018393
At iterate    92  f =       80.936  |proj g|=      0.020207
At iterate    93  f =       80.935  |proj g|=      0.050297
At iterate    94  f =       80.935  |proj g|=      0.018193
At iterate    95  f =       80.935  |proj g|=      0.015472
At iterate    96  f =       80.935  |proj g|=      0.014675
At iterate    97  f =       80.935  |proj g|=      0.072348
At iterate    98  f =       80.935  |proj g|=       0.02166
At iterate    99  f =       80.935  |proj g|=      0.017904
At iterate   100  f =       80.935  |proj g|=      0.040947
At iterate   101  f =       80.935  |proj g|=      0.024479
final  value 80.934858 
stopped after 101 iterations

5.4 Augment the design.

The function addNroyDesignPoints builds an emulator for each model output in Y. It compares the output of each emulator at a number of candidate desin points, and chooses a space-filling set of them that that are Not Ruled Out Yet (statistically close to the observation at Y_target).

5.5 Write the augmented design

The function write_jules_design here simply takes the points calculated by addNroyDesignPoints and writes them to configuration files.

5.6 Check the design

Check that the augmented design produces what we expect. New ensemble members should be somewhat constrained within the boundaries of the original design, if the comparison to data offers any constraint.

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -109.2909 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       109.29  |proj g|=        1.871
At iterate     1  f =       100.41  |proj g|=        1.3975
At iterate     2  f =       96.533  |proj g|=        1.3088
At iterate     3  f =       95.227  |proj g|=        1.2724
At iterate     4  f =       92.512  |proj g|=        1.5964
At iterate     5  f =       91.143  |proj g|=        1.7368
At iterate     6  f =       89.279  |proj g|=        1.6991
At iterate     7  f =       88.195  |proj g|=        1.7046
At iterate     8  f =       87.014  |proj g|=        1.5114
At iterate     9  f =       85.689  |proj g|=        1.0118
At iterate    10  f =       85.023  |proj g|=        1.7091
At iterate    11  f =       84.569  |proj g|=        1.3563
At iterate    12  f =       84.158  |proj g|=       0.80825
At iterate    13  f =        84.09  |proj g|=        1.6904
At iterate    14  f =       83.823  |proj g|=        1.6909
At iterate    15  f =        83.55  |proj g|=       0.81435
At iterate    16  f =       83.371  |proj g|=       0.96783
At iterate    17  f =       83.094  |proj g|=        1.1149
At iterate    18  f =       82.902  |proj g|=       0.63225
At iterate    19  f =       82.876  |proj g|=        1.6631
At iterate    20  f =       82.821  |proj g|=       0.47627
At iterate    21  f =       82.812  |proj g|=       0.44135
At iterate    22  f =       82.724  |proj g|=       0.87459
At iterate    23  f =       82.656  |proj g|=       0.80084
At iterate    24  f =       82.618  |proj g|=       0.37465
At iterate    25  f =       82.608  |proj g|=       0.27736
At iterate    26  f =       82.599  |proj g|=       0.18904
At iterate    27  f =       82.587  |proj g|=       0.29552
At iterate    28  f =       82.563  |proj g|=       0.53817
At iterate    29  f =       82.536  |proj g|=         0.389
At iterate    30  f =        82.52  |proj g|=       0.27542
At iterate    31  f =       82.509  |proj g|=       0.43571
At iterate    32  f =       82.496  |proj g|=       0.39767
At iterate    33  f =       82.472  |proj g|=       0.67351
At iterate    34  f =       82.447  |proj g|=        1.0541
At iterate    35  f =       82.426  |proj g|=       0.41252
At iterate    36  f =       82.407  |proj g|=       0.36104
At iterate    37  f =        82.39  |proj g|=       0.37105
At iterate    38  f =       82.352  |proj g|=       0.40153
At iterate    39  f =       82.345  |proj g|=       0.48576
At iterate    40  f =       82.338  |proj g|=       0.30354
At iterate    41  f =       82.333  |proj g|=       0.28668
At iterate    42  f =       82.321  |proj g|=       0.24732
At iterate    43  f =       82.309  |proj g|=       0.34561
At iterate    44  f =       82.295  |proj g|=       0.29314
At iterate    45  f =       82.283  |proj g|=        0.5463
At iterate    46  f =       82.269  |proj g|=       0.55403
At iterate    47  f =       82.266  |proj g|=       0.38142
At iterate    48  f =       82.264  |proj g|=       0.15815
At iterate    49  f =        82.26  |proj g|=       0.11062
At iterate    50  f =       82.258  |proj g|=        0.1582
At iterate    51  f =       82.255  |proj g|=       0.22033
At iterate    52  f =       82.254  |proj g|=       0.13105
At iterate    53  f =       82.253  |proj g|=       0.05804
At iterate    54  f =       82.253  |proj g|=       0.15167
At iterate    55  f =       82.253  |proj g|=      0.054998
At iterate    56  f =       82.252  |proj g|=       0.03075
At iterate    57  f =       82.252  |proj g|=      0.029696
At iterate    58  f =       82.252  |proj g|=      0.016543
At iterate    59  f =       82.252  |proj g|=      0.014205
At iterate    60  f =       82.252  |proj g|=      0.026848
At iterate    61  f =       82.252  |proj g|=      0.027097
At iterate    62  f =       82.252  |proj g|=      0.012538
At iterate    63  f =       82.252  |proj g|=      0.011962
At iterate    64  f =       82.252  |proj g|=      0.011229
At iterate    65  f =       82.252  |proj g|=       0.01993
At iterate    66  f =       82.252  |proj g|=      0.012827
At iterate    67  f =       82.252  |proj g|=      0.016769
At iterate    68  f =       82.252  |proj g|=      0.038307
At iterate    69  f =       82.252  |proj g|=       0.02021
At iterate    70  f =       82.252  |proj g|=       0.01395
At iterate    71  f =       82.252  |proj g|=      0.019049
At iterate    72  f =       82.252  |proj g|=      0.028594
At iterate    73  f =       82.252  |proj g|=      0.033337
At iterate    74  f =       82.252  |proj g|=      0.020535
At iterate    75  f =       82.252  |proj g|=     0.0075909
At iterate    76  f =       82.252  |proj g|=     0.0055178
At iterate    77  f =       82.252  |proj g|=      0.033389
At iterate    78  f =       82.252  |proj g|=     0.0071442
At iterate    79  f =       82.252  |proj g|=        0.0111
At iterate    80  f =       82.252  |proj g|=       0.01729
At iterate    81  f =       82.252  |proj g|=       0.01571
At iterate    82  f =       82.252  |proj g|=      0.029493
At iterate    83  f =       82.252  |proj g|=     0.0070285
At iterate    84  f =       82.252  |proj g|=     0.0084202
At iterate    85  f =       82.252  |proj g|=      0.012163
At iterate    86  f =       82.252  |proj g|=      0.089832
At iterate    87  f =       82.252  |proj g|=       0.03557
At iterate    88  f =       82.252  |proj g|=      0.010387
At iterate    89  f =       82.252  |proj g|=       0.02101
At iterate    90  f =       82.252  |proj g|=      0.017472
At iterate    91  f =       82.252  |proj g|=      0.025682
At iterate    92  f =       82.252  |proj g|=     0.0079639
At iterate    93  f =       82.252  |proj g|=     0.0073904
At iterate    94  f =       82.252  |proj g|=      0.018008
At iterate    95  f =       82.252  |proj g|=     0.0062852
At iterate    96  f =       82.252  |proj g|=     0.0046032
At iterate    97  f =       82.252  |proj g|=     0.0028965
At iterate    98  f =       82.251  |proj g|=     0.0083294
At iterate    99  f =       82.251  |proj g|=      0.004218
At iterate   100  f =       82.251  |proj g|=     0.0024721
At iterate   101  f =       82.251  |proj g|=     0.0026606
final  value 82.251495 
stopped after 101 iterations

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -1383.375 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1383.4  |proj g|=       1.7664
At iterate     1  f =       1372.9  |proj g|=        1.6648
At iterate     2  f =       1364.8  |proj g|=        1.7747
At iterate     3  f =       1356.7  |proj g|=        1.8961
At iterate     4  f =       1350.1  |proj g|=        1.7144
At iterate     5  f =       1347.8  |proj g|=        1.6582
At iterate     6  f =       1341.7  |proj g|=         1.481
At iterate     7  f =       1337.6  |proj g|=        1.6235
At iterate     8  f =       1335.4  |proj g|=        1.4778
At iterate     9  f =         1334  |proj g|=       0.92855
At iterate    10  f =         1332  |proj g|=        1.5095
At iterate    11  f =       1330.6  |proj g|=        1.0917
At iterate    12  f =       1329.3  |proj g|=        1.2354
At iterate    13  f =       1328.5  |proj g|=       0.82781
At iterate    14  f =       1328.3  |proj g|=        0.9746
At iterate    15  f =       1328.1  |proj g|=       0.55288
At iterate    16  f =         1328  |proj g|=       0.39538
At iterate    17  f =       1327.8  |proj g|=       0.30206
At iterate    18  f =       1327.8  |proj g|=        1.4645
At iterate    19  f =       1327.8  |proj g|=       0.56521
At iterate    20  f =       1327.8  |proj g|=       0.21692
At iterate    21  f =       1327.8  |proj g|=       0.56111
At iterate    22  f =       1327.8  |proj g|=       0.51424
At iterate    23  f =       1327.8  |proj g|=       0.53842
At iterate    24  f =       1327.8  |proj g|=        0.1763
At iterate    25  f =       1327.7  |proj g|=       0.20315
At iterate    26  f =       1327.7  |proj g|=       0.14516
At iterate    27  f =       1327.7  |proj g|=      0.054732
At iterate    28  f =       1327.7  |proj g|=      0.084036
At iterate    29  f =       1327.7  |proj g|=      0.045267
At iterate    30  f =       1327.7  |proj g|=      0.046735
At iterate    31  f =       1327.7  |proj g|=       0.14927
At iterate    32  f =       1327.7  |proj g|=       0.11387
At iterate    33  f =       1327.7  |proj g|=      0.025094
At iterate    34  f =       1327.7  |proj g|=      0.018372
At iterate    35  f =       1327.7  |proj g|=      0.018539
At iterate    36  f =       1327.7  |proj g|=      0.017334
At iterate    37  f =       1327.7  |proj g|=      0.035245
At iterate    38  f =       1327.7  |proj g|=      0.016258
At iterate    39  f =       1327.7  |proj g|=     0.0075892
At iterate    40  f =       1327.7  |proj g|=      0.008652
At iterate    41  f =       1327.7  |proj g|=      0.023457

iterations 41
function evaluations 48
segments explored during Cauchy searches 62
BFGS updates skipped 0
active bounds at final generalized Cauchy point 20
norm of the final projected gradient 0.0234567
final function value 1327.74

F = 1327.74
final  value 1327.742870 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -2561.221 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       2561.2  |proj g|=       1.8661
At iterate     1  f =       2533.6  |proj g|=             2
At iterate     2  f =       2519.6  |proj g|=             2
At iterate     3  f =       2516.5  |proj g|=        1.9534
At iterate     4  f =       2510.4  |proj g|=        1.7867
At iterate     5  f =         2505  |proj g|=        1.6826
At iterate     6  f =       2501.5  |proj g|=         1.618
At iterate     7  f =       2494.2  |proj g|=        1.3844
At iterate     8  f =       2492.2  |proj g|=        1.6032
At iterate     9  f =       2491.3  |proj g|=         1.589
At iterate    10  f =       2488.9  |proj g|=          1.53
At iterate    11  f =       2487.5  |proj g|=        1.0728
At iterate    12  f =       2487.2  |proj g|=        1.0711
At iterate    13  f =         2487  |proj g|=        1.1027
At iterate    14  f =       2486.8  |proj g|=       0.47609
At iterate    15  f =       2486.8  |proj g|=       0.42024
At iterate    16  f =       2486.8  |proj g|=        1.2305
At iterate    17  f =       2486.8  |proj g|=       0.38381
At iterate    18  f =       2486.8  |proj g|=       0.22702
At iterate    19  f =       2486.8  |proj g|=       0.21667
At iterate    20  f =       2486.8  |proj g|=       0.17646
At iterate    21  f =       2486.7  |proj g|=         0.202
At iterate    22  f =       2486.7  |proj g|=       0.44121
At iterate    23  f =       2486.7  |proj g|=       0.16229
At iterate    24  f =       2486.7  |proj g|=      0.066349
At iterate    25  f =       2486.7  |proj g|=      0.018863
At iterate    26  f =       2486.7  |proj g|=      0.034152
At iterate    27  f =       2486.7  |proj g|=      0.043297
At iterate    28  f =       2486.7  |proj g|=      0.014902
At iterate    29  f =       2486.7  |proj g|=      0.016268
At iterate    30  f =       2486.7  |proj g|=      0.013669
At iterate    31  f =       2486.7  |proj g|=      0.021967

iterations 31
function evaluations 34
segments explored during Cauchy searches 54
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0219674
final function value 2486.74

F = 2486.74
final  value 2486.744525 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io + 
    dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io + 
    g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io + 
    lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io + 
    retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io + 
    sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -2315.205 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       2315.2  |proj g|=       1.9221
At iterate     1  f =       2314.2  |proj g|=        1.7645
At iterate     2  f =       2314.1  |proj g|=       0.74274
At iterate     3  f =       2313.9  |proj g|=        1.3088
At iterate     4  f =       2313.7  |proj g|=        1.9298
At iterate     5  f =       2313.5  |proj g|=        1.9378
At iterate     6  f =       2313.1  |proj g|=        1.9199
At iterate     7  f =         2312  |proj g|=        1.9412
At iterate     8  f =       2311.8  |proj g|=        1.9461
At iterate     9  f =       2310.4  |proj g|=        1.3481
At iterate    10  f =       2309.3  |proj g|=        1.2924
At iterate    11  f =       2308.7  |proj g|=         1.937
At iterate    12  f =       2308.3  |proj g|=        1.5008
At iterate    13  f =       2308.3  |proj g|=         1.216
At iterate    14  f =       2308.2  |proj g|=        1.8896
At iterate    15  f =       2308.1  |proj g|=        1.8935
At iterate    16  f =       2307.5  |proj g|=        1.9032
At iterate    17  f =         2307  |proj g|=        1.8988
At iterate    18  f =       2304.7  |proj g|=        1.0033
At iterate    19  f =       2303.6  |proj g|=        1.5714
At iterate    20  f =       2303.5  |proj g|=         1.917
At iterate    21  f =       2303.3  |proj g|=        1.9143
At iterate    22  f =       2303.2  |proj g|=        1.9086
At iterate    23  f =       2303.1  |proj g|=        1.8843
At iterate    24  f =       2302.9  |proj g|=       0.70255
At iterate    25  f =       2302.9  |proj g|=       0.68074
At iterate    26  f =       2302.9  |proj g|=        1.9109
At iterate    27  f =       2302.8  |proj g|=        1.7663
At iterate    28  f =       2302.8  |proj g|=       0.67488
At iterate    29  f =       2302.6  |proj g|=        1.0907
At iterate    30  f =       2302.3  |proj g|=        1.9095
At iterate    31  f =       2302.2  |proj g|=        1.9105
At iterate    32  f =         2302  |proj g|=        1.4736
At iterate    33  f =         2302  |proj g|=        1.4695
At iterate    34  f =       2301.8  |proj g|=          1.32
At iterate    35  f =       2301.7  |proj g|=          1.91
At iterate    36  f =       2301.7  |proj g|=        1.3149
At iterate    37  f =       2301.6  |proj g|=       0.83873
At iterate    38  f =       2301.5  |proj g|=       0.96389
At iterate    39  f =       2301.5  |proj g|=        1.3464
At iterate    40  f =       2301.4  |proj g|=       0.53138
At iterate    41  f =       2301.3  |proj g|=        1.8834
At iterate    42  f =       2301.3  |proj g|=       0.82272
At iterate    43  f =       2301.3  |proj g|=       0.75558
At iterate    44  f =       2301.3  |proj g|=        1.8803
At iterate    45  f =       2301.3  |proj g|=        1.8814
At iterate    46  f =       2301.3  |proj g|=        1.8812
At iterate    47  f =       2301.3  |proj g|=         1.881
At iterate    48  f =       2301.2  |proj g|=        1.3784
At iterate    49  f =       2301.2  |proj g|=       0.62431
At iterate    50  f =       2301.2  |proj g|=        1.9114
At iterate    51  f =       2301.2  |proj g|=        1.9123
At iterate    52  f =       2301.1  |proj g|=        1.9116
At iterate    53  f =       2301.1  |proj g|=       0.63241
At iterate    54  f =       2301.1  |proj g|=       0.28367
At iterate    55  f =       2301.1  |proj g|=       0.31082
At iterate    56  f =       2301.1  |proj g|=       0.62345
At iterate    57  f =       2301.1  |proj g|=         1.911
At iterate    58  f =       2301.1  |proj g|=        1.9119
At iterate    59  f =       2301.1  |proj g|=        1.9112
At iterate    60  f =         2301  |proj g|=       0.29321
At iterate    61  f =         2301  |proj g|=       0.31545
At iterate    62  f =         2301  |proj g|=       0.84934
At iterate    63  f =         2301  |proj g|=        1.8836
At iterate    64  f =         2301  |proj g|=        1.8821
At iterate    65  f =         2301  |proj g|=       0.58887
At iterate    66  f =         2301  |proj g|=        1.2003
At iterate    67  f =         2301  |proj g|=       0.43465
At iterate    68  f =         2301  |proj g|=       0.42034
At iterate    69  f =         2301  |proj g|=       0.60752
At iterate    70  f =         2301  |proj g|=        1.6659
At iterate    71  f =         2301  |proj g|=        1.8816
At iterate    72  f =       2300.9  |proj g|=        1.8816
At iterate    73  f =       2300.9  |proj g|=        1.8828
At iterate    74  f =       2300.9  |proj g|=        1.8811
At iterate    75  f =       2300.9  |proj g|=       0.40944
At iterate    76  f =       2300.9  |proj g|=       0.36834
At iterate    77  f =       2300.9  |proj g|=       0.74183
At iterate    78  f =       2300.9  |proj g|=        1.6322
At iterate    79  f =       2300.9  |proj g|=         1.914
At iterate    80  f =       2300.9  |proj g|=        1.9141
At iterate    81  f =       2300.9  |proj g|=        1.9155
At iterate    82  f =       2300.9  |proj g|=        1.1075
At iterate    83  f =       2300.8  |proj g|=       0.97412
At iterate    84  f =       2300.8  |proj g|=       0.96614
At iterate    85  f =       2300.8  |proj g|=        1.8816
At iterate    86  f =       2300.8  |proj g|=        1.8815
At iterate    87  f =       2300.8  |proj g|=        1.9143
At iterate    88  f =       2300.7  |proj g|=       0.61654
At iterate    89  f =       2300.7  |proj g|=       0.61993
At iterate    90  f =       2300.7  |proj g|=       0.61744
At iterate    91  f =       2300.7  |proj g|=       0.59652
At iterate    92  f =       2300.7  |proj g|=       0.47719
At iterate    93  f =       2300.7  |proj g|=       0.53494
At iterate    94  f =       2300.7  |proj g|=       0.53394
At iterate    95  f =       2300.7  |proj g|=       0.52042
At iterate    96  f =       2300.7  |proj g|=       0.60048
At iterate    97  f =       2300.7  |proj g|=        1.9142
At iterate    98  f =       2300.7  |proj g|=        1.9152
At iterate    99  f =       2300.6  |proj g|=        1.9161
At iterate   100  f =       2300.6  |proj g|=        1.9168
At iterate   101  f =       2300.6  |proj g|=        1.3656
final  value 2300.595796 
stopped after 101 iterations

Visualise the predicted outputs at the NROY points of the old design, and the new suggested design.

5.7 Visualising the emulated outputs at the proposed new design points.

Interestingly, these aren’t perfectly within the original hard boundaries set by Andy, even though I’ve set those boundaries to be the +- 4 standard deviation threholds in the History Match. I suggest this is because there is model discrepancy, and that there is considerable wriggle room induced from emulator uncertainty.

In particular, it appears that vegetation carbon is difficult to keep high, and that many NROY proposed members have a fairly low vegetation carbon. This might need a discrepancy term, or adjusting in some other way. It certainly needs exploring, and a OAAT plot might give clues as to the parameters to choose.

How good are the four emulators that we’ve built? Are there biases? (there’s no real evidence of this)

5.8 One-at-a-time sensitivity analysis for understanding model response

It’s obviously hard to maintain a high vegetation carbon in particular. What parameter values might you choose to do this, and what might be the trade-offs you have to make?

First, what parameters affect vegetation carbon and how? How sure are we about that?

some values will be clipped

What might be the trade-offs for a high (or accurate) vegetation carbon? are they acceptable? Plot the oaat sensitivity of the other 3 outputs we’re calibrating on.

What do the emulators make of the design points which actually make Andy’s “hard boundary” criteria? If we leave them out, do they still place the output within the hard boundaries?

5.9 Check the emulators that produce the new design

Do a leave-one-out cross validation of points inside the hard boundaries using the wave1 fits.

We see in the leave-one-out analysis that the emulator is consistently under-predicting the vegetation carbon (though the uncertainty estimate often covers the actual value).This suggests (1) that there isn’t really a huge problem with a model discrepancy (or at least that isn’t the only problem), and (2) the history matching is working as it should, and taking into account a not-great emulator.

5.10 Compare the straight km with a two-step emulator

Is a two-step emulator any better at emulating those crucial points which fall within aw’s hard boundaries? First, create a list of emulator fits.

It doesn’t appear that the two-step emulator (here plotted in red) is doing any better than the regular emulator.

Two other things I can think of to check: 1) how about using “multistart” to choose different starting conditions for optimising the emulators and 2) Using a flat prior for the mean function rather than a linear prior.

Loading required package: Matrix

Attaching package: ‘Matrix’

The following object is masked from ‘package:spam’:

    det

Loaded glmnet 4.1-1

6 Exploring new constraints

Input from Eddy Robertson:

“We can assume that the majority of vegetation carbon is stored in tree trunks so, the carbon density of trees is approximately (cVeg_lnd_mean / treeFrac_lnd_mean). Although I suppose it’s possible that in the PPE the shrubs have become tree-like, so it might be interesting to plot cVeg_lnd_mean versus (shrubFrac_lnd_mean + treeFrac_lnd_mean) as well. The first question is whether cVeg_lnd_mean increases linearly with treeFrac_lnd_mean. Then if this turned out to be interesting, we could produce a treeCVeg_lnd_mean output.”

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~1
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -324.7357 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       324.74  |proj g|=       1.8894
At iterate     1  f =       154.68  |proj g|=             2
At iterate     2  f =       144.37  |proj g|=        1.9477
At iterate     3  f =       140.12  |proj g|=        1.8272
At iterate     4  f =       137.65  |proj g|=        1.5752
At iterate     5  f =       135.35  |proj g|=        1.4243
At iterate     6  f =       134.81  |proj g|=        1.1667
At iterate     7  f =       134.52  |proj g|=        1.0068
At iterate     8  f =       134.47  |proj g|=       0.55475
At iterate     9  f =       134.44  |proj g|=       0.37207
At iterate    10  f =       134.43  |proj g|=       0.20987
At iterate    11  f =       134.42  |proj g|=      0.076796
At iterate    12  f =       134.42  |proj g|=      0.026753
At iterate    13  f =       134.42  |proj g|=       0.01759
At iterate    14  f =       134.42  |proj g|=     0.0060535
At iterate    15  f =       134.42  |proj g|=     0.0041787
At iterate    16  f =       134.42  |proj g|=     0.0021465
At iterate    17  f =       134.42  |proj g|=       0.00377
At iterate    18  f =       134.42  |proj g|=      0.001123
At iterate    19  f =       134.42  |proj g|=    0.00021622

iterations 19
function evaluations 20
segments explored during Cauchy searches 42
BFGS updates skipped 0
active bounds at final generalized Cauchy point 19
norm of the final projected gradient 0.000216221
final function value 134.421

F = 134.421
final  value 134.420708 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~1
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -1621.762 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       1621.8  |proj g|=       1.7706
At iterate     1  f =       1507.2  |proj g|=        1.4766
At iterate     2  f =       1423.3  |proj g|=        1.7768
At iterate     3  f =       1418.6  |proj g|=        1.4019
At iterate     4  f =       1415.9  |proj g|=        1.6446
At iterate     5  f =       1415.3  |proj g|=        1.7139
At iterate     6  f =       1415.1  |proj g|=        1.6554
At iterate     7  f =       1414.8  |proj g|=        1.4108
At iterate     8  f =       1414.7  |proj g|=       0.99795
At iterate     9  f =       1414.5  |proj g|=       0.73402
At iterate    10  f =       1414.5  |proj g|=       0.66567
At iterate    11  f =       1414.5  |proj g|=       0.13919
At iterate    12  f =       1414.5  |proj g|=        0.0858
At iterate    13  f =       1414.5  |proj g|=      0.034268
At iterate    14  f =       1414.5  |proj g|=      0.063708
At iterate    15  f =       1414.5  |proj g|=      0.010821
At iterate    16  f =       1414.5  |proj g|=     0.0075671
At iterate    17  f =       1414.5  |proj g|=      0.012372
At iterate    18  f =       1414.5  |proj g|=     0.0071278
At iterate    19  f =       1414.5  |proj g|=     0.0036723

iterations 19
function evaluations 21
segments explored during Cauchy searches 70
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.00367225
final function value 1414.5

F = 1414.5
final  value 1414.500017 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~1
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -2800.645 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=       2800.6  |proj g|=       1.8463
At iterate     1  f =       2616.7  |proj g|=        1.6218
At iterate     2  f =       2605.4  |proj g|=         1.553
At iterate     3  f =       2583.7  |proj g|=        1.6224
At iterate     4  f =       2577.2  |proj g|=         1.011
At iterate     5  f =         2573  |proj g|=        1.3781
At iterate     6  f =       2571.3  |proj g|=         1.308
At iterate     7  f =       2571.1  |proj g|=        0.9583
At iterate     8  f =         2571  |proj g|=       0.62272
At iterate     9  f =       2570.9  |proj g|=       0.44321
At iterate    10  f =       2570.9  |proj g|=       0.12021
At iterate    11  f =       2570.9  |proj g|=      0.055708
At iterate    12  f =       2570.9  |proj g|=      0.016351
At iterate    13  f =       2570.9  |proj g|=     0.0027779

iterations 13
function evaluations 16
segments explored during Cauchy searches 52
BFGS updates skipped 0
active bounds at final generalized Cauchy point 27
norm of the final projected gradient 0.00277786
final function value 2570.93

F = 2570.93
final  value 2570.925355 
converged

optimisation start
------------------
* estimation method   : MLE 
* optimisation method : BFGS 
* analytical gradient : used
* trend model : ~1
* covariance model : 
  - type :  matern5_2 
  - nugget : NO
  - parameters lower bounds :  1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 
  - parameters upper bounds :  2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2 
  - best initial criterion value(s) :  -2465.963 

N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate     0  f=         2466  |proj g|=       1.7876
At iterate     1  f =       2357.6  |proj g|=             2
At iterate     2  f =       2342.2  |proj g|=             2
At iterate     3  f =       2334.4  |proj g|=        1.9202
At iterate     4  f =       2331.2  |proj g|=        1.8804
At iterate     5  f =         2325  |proj g|=        1.7272
At iterate     6  f =         2322  |proj g|=        1.6469
At iterate     7  f =       2320.5  |proj g|=        1.4532
At iterate     8  f =       2317.5  |proj g|=        1.4216
At iterate     9  f =         2317  |proj g|=        1.3362
At iterate    10  f =       2316.1  |proj g|=        1.2141
At iterate    11  f =       2315.5  |proj g|=        1.1842
At iterate    12  f =       2315.4  |proj g|=       0.93199
At iterate    13  f =       2315.4  |proj g|=        0.5059
At iterate    14  f =       2315.4  |proj g|=       0.38256
At iterate    15  f =       2315.4  |proj g|=       0.53325
At iterate    16  f =       2315.4  |proj g|=       0.18709
At iterate    17  f =       2315.4  |proj g|=      0.038636
At iterate    18  f =       2315.4  |proj g|=       0.02819
At iterate    19  f =       2315.4  |proj g|=      0.033483
At iterate    20  f =       2315.4  |proj g|=      0.015972
At iterate    21  f =       2315.4  |proj g|=      0.014807
At iterate    22  f =       2315.4  |proj g|=     0.0076813
At iterate    23  f =       2315.4  |proj g|=     0.0071707

iterations 23
function evaluations 27
segments explored during Cauchy searches 50
BFGS updates skipped 0
active bounds at final generalized Cauchy point 21
norm of the final projected gradient 0.00717074
final function value 2315.38

F = 2315.38
final  value 2315.375074 
converged

Is the range of carbon densities narrower in the “hard boundary” set?

6.1 Now that we have constraints, we can explore their effect on the output

6.2 How much does constraining an output constrain other outputs?

For example, what effect has the AW “hard boundary” constraint produced on ranges of the other (marginal) outputs? Which constraint is doing the majority of the work?

A function to find the proportion of output space that is removed (or retained) when applying a constraint.

All constraints together take the normalized range of all the outputs to a little under a half of the original range, on average. Of course, four of the inputs are constrained here, which has quite a large impact.

We can see what impacts the constraints have individually

The NBP constraint has no marginal impact

Vegetation carbon (cVeg) has the strongest impact on everything else. This is no suprise, given that so few of the ensemble members seem to simulate carbon vegetation high enough.

6.3 Writing robust, parallelised, multistart emulators

7 Direct twoStep emulator using foreach